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Abstract. We present a local algorithm (constant-time distributed algorithm) for ap- 
proximating max-min LPs. The objective is to maximise u> subject to Ax < 1, Cx > loI, 
and X > for nonnegative matrices A and C. The approximation ratio of our algorithm is 
the best possible for any local algorithm; there is a matching unconditional lower bound. 



1. Introduction 

In a max-min linear program (max-min LP), the objective is to 

maximise min^g^^ c^x 

subject to Ax < 1, (1.1) 

X > 

or, equivalently, to maximise uj subject to Ax < 1, Cx > ujI, and x > 0. The matrices A 
and C are nonnegative and sparse: each row of A has at most A/ positive elements, and 
each row of C has at most Ak positive elements; here A/ and A^^- are constants. 

Our work studies local algorithms [HIIIO], that is, distributed algorithms that complete 
in constant time (constant number of synchronous communication rounds), independent 
of the size of the network. We assume that there is a network with one node v € V for 
each variable Xy, one node i I for each constraint a^x < 1, and one node k (z K for 
each objective c^x. Nodes v (z V and i (z I are adjacent if is positive, and nodes 
V (zV and k ^ K are adjacent if Ckv is positive. In a local algorithm for max-min LPs, the 
node V must choose the value based on the information that is available within its 
constant-radius neighbourhood in the network. 

Max-min LPs are a generalisation of packing LPs. Direct applications of max-min 
LPs include various tasks of fair resource allocation, such as fair bandwidth allocation in 
a communication network and balanced data gathering in a wireless sensor network. An 
algorithm for approximating max-min LPs allows one to solve approximate mixed packing 
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and covering problems |12j as well. Special cases of mixed packing and covering LPs include 
finding an (approximate) solution to a nonnegative system of linear equations. 

1.1. Prior work 

For any e > 0, there exists a local (1 + e)-approximation algorithm for packing and 
covering LPs, assuming a bounded-degree graph and bounded coefficients [318]. However, 
this is not the case with max-min LPs: 

Theorem 1.1 (Floreen et al. [4]). For any A/ > 2 and Ak > 2, there exists no local approx- 
imation algorithm for the max-min LP problem with the approximation ratio A/(l — 1/Aj^). 

This theorem holds even in the following special cases: (i) bipartite max-min LPs, where 
each column of A and each column of C contains only one nonzero element, and (ii) 0/1 
max-min LPs, where each element of A and C is either or 1. 

The approximation ratio A/(l — l/A^) + e for any positive e can be achieved by a local 
algorithm for bipartite problems [1] and for a class of 0/1 problems [3]. However, for general 
max-min LPs, the best known local algorithm has been the safe algorithm that achieves 
the factor A/ approximation [Sllllj. 

1.2. Contribution 

The following theorem summarises the main contribution of this paper. 

Theorem 1.2. For any A/ > 2, Ak > 2, and e > 0, there exists a local approximation 
algorithm for the max-min LP problem with the approximation ratio A/(l — 1/ Ak) + £• 

By Theorem II. H the algorithm in Theorem 11.21 is optimal in the sense that no local 
algorithm can achieve a better approximation ratio. For the case Aj = Ak = 2, the 
approximation ratio is 1 -|- e, that is, there is a local approximation scheme. 

1.3. Definitions 

We now give a detailed definition of a max-min LP in a distributed setting. Let Q = 
(yuLUK,E) be a bipartite, undirected communication graph. The nodes v £ V are called 
agents, the nodes i (z L are called constraints, and the nodes k £ K are called objectives; 
the sets V , L, and K are disjoint. Edges e G E are of the form e = {v, i} or e = {v, k} 
where v & V , i (z L , and k (z K. 

Let Vi = {v eV : {v, i} e E}, Vk = {v e V : {v, k} £ E}, L^ = {i £ L : {v, i} £ E}, and 

= {k £ K : {v, k} £ E} for alii £ L, k £ K, and v £V. We assume that \Vi\ < Aj and 

< Ak for alH G / and k £ K for some constants A/ and Ak- 

A max-min linear program associated with Q is defined as follows. Associate a variable 
Xy with each agent v £V , associate a coefficient a^^ > with each edge {i, v} £ E, i £ L, 
V £ V, and associate a coefficient c^v > with each edge {k, v}£E,k£K,v£V. The 
task is to 

maximise uj{x) = min^gA' ^veVk ^^vXv 

subject to Yliv&v^ (^iv^v < 1, yi£ L, (1.2) 

Xv > 0, Vt; e V. 

The local input of an agent v £ V consists of the sets and and the coefficients 
o-iv^Ckv for aW. i £ Ly,k £ Ky. The local input of a constraint i £ L consists of Vi, and the 
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local input of an objective k £ K consists of 14 . In a local algorithm, there is a constant D 
such that each agent v G V chooses the output Xy based on the local inputs of the nodes 
within distance D (in number of edges) from v in the communication graph Q; the constant 
D is the local horizon of the algorithm. 

Theorem 11.11 holds even if we assume that each node of the graph Q has a unique 
identifiers. We show that Theorem 11.21 holds even if we do not have unique identifiers. We 
merely assume port numbering [Ij: each edge {s, t} in Q has two natural numbers associated 
with it, the port number in s and the port number in t. 

1.4. Overview of the algorithm 

We begin in f|2]by reminding that, in the port numbering model, we can without loss 
of generality focus on the case where the communication graph Q is a (countably infinite 
but locally finite) tree [HHj. 

In f|3l we present a series of local transformations that simplify the structure of the 
problem. We show that with the objective of establishing Theorem \1.2\ it is sufficient to 
focus on the special case where \Vi\ = 2, \Vk\ > 2, \Ky\ = 1, > 1, and c^u = 1 for all 
i (z I , k & K , V (z V , u (z Vk- Figure [T] shows an example of a communication graph G after 
the local transformations. 

In ^ we present a local algorithm for this special case. In ^ we prove that the output 
of the algorithm is a factor 2(1 — 1/Ak) + e' approximation. Put together, we have a local, 
factor A/(l — 1/Ak) + e approximation for general max-min LPs, for any e > 0. 

The intuition behind the algorithm in ^ is best understood if we study its analysis 
in ^ In the analysis, it is convenient to assume that we have assigned a one-dimensional 
coordinate, layer, to each node of the tree G; see Figure [J for an example. When we 
assign the layers, we also partition the agents into up- agents and down- agents. We have 
alternatingly layers of up-agents, constraints, down-agents, and objectives. Each objective 
has exactly one adjacent up-agent "above" it, and at least one adjacent down-agent "below" 
it; hence the names "up" and "down". 

We now assign every R^^ layer of objectives to be passive, including the adjacent agents 
that set Xy = 0. Each agent v € V computes an upper bound t^ of the optimum; see §4.21 
Then we construct a solution for the active layers in a greedy manner, starting with a layer 
of passive up-agents and propagating information upwards until we reach the next layer of 
passive down-agents; see §4.31 

This is not yet an approximation for the original problem: while most objectives perform 
at least as well as in the global optimum, the passive objectives have utility 0. By applying 
ideas from the shifting strategy [2l[6], we could consider R possible choices for the locations 
of the passive layers. Then we could take averages over these to obtain a solution y. In 
§5.21 we show that y would indeed be a factor R/{R — 1) approximation. 

There is one difficulty, however. We cannot assign the layers by a local algorithm in a 
globally consistent manner; in particular, we do not know whether a given agent is a down- 
agent or an up-agent. To overcome this, we consider both possible roles for each agent, up 
and down. For both roles, we compute a candidate solution by applying the shifting strategy. 
Finally we take the average of both candidate solutions. This is the essence of (|4.13p . 

In ^ we prove that this local approach yields a globally feasible solution, and the 
solution is within factor 2(1 — 1/Aj^) -|- e' of the optimum. The constant e' > can be 
made arbitrarily small by choosing a sufficiently large R. 
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Figure 1: The graph Q (after the local transformations in ^ and the layers (see We 
have chosen R = 3 and hence r = 1. Thick lines highlight the tree Au (see y4.ip . 
There are several ways to choose the layers; nevertheless, if u is an up-agent on 
layer —1 then the levels in Au necessarily coincide with the layers. 



Yet another basic hurdle implicitly overcome in the proof of Theorem 11.21 stems from 
underconstrained instances: if there are several equally good solutions, one needs to choose 
between them in a globally consistent manner. Our definition of the values g'^ and g~ 
in ()4.8p - (l4.10p addresses this by focusing on a particular extreme point. Each layer of 
down-agents chooses as large values as possible, without violating the constraints "be- 
low" them. Each layer of up-agents chooses as small values g~ as possible, as long as the 
objectives "below" them meet the smoothed upper bounds Sv 



2. Port numbering and unfolding 

Our algorithm does not need to use any node identifiers; port numbering is sufficient. 
In the port numbering model, a local algorithm cannot distinguish between a short cycle 
and an infinitely long path. We can exploit this limitation to simplify the description of our 
local algorithm: we can assume that we have unfolded all cycles of the graph Q y4|4j. 

The definition of unfolding requires some preliminaries. A walk of length i in a graph 
^ is a nonempty tuple (uq, ei, mi, 62, M2, • • • , e^, ui) of alternating nodes and edges in G such 
that, for all j = 1,2, . . . , I, the edge ej joins the nodes Uj-i and uj. The walk is said to start 
at uq and end at n^. A walk is non-backtracking if ej-i 7^ ej holds for all j = 1, 2, . . . A 
path is a walk with no repeated nodes. 

Let ^ be a finite connected graph and let r be a node of Q. The unfolding of Q rooted 
at r is the undirected simple graph Q' obtained as follows. The node set of Q' is the set of 
all non-backtracking walks in Q that start at r. Two nodes of G' are joined by an edge iff 
one can be obtained from the other by appending exactly one edge and one node of Q. 
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We associate with each node of Q' a parent node of namely the end-node of the walk. 
We also associate with each edge of Q' a parent edge of t/, namely the appended edge. 

Remarks. 

(1) The unfolding Q' is a tree. 

(2) The unfolding Q' is finite iff ^ is a tree; otherwise Q' is countably infinite. 

(3) Any two unfoldings of Q rooted at different nodes are isomorphic. In what follows 
we refer to "the" unfolding of Q without specifying a particular root node. 

(4) Assuming that the graph Q has port numbers associated with the ends of its edges, 
the unfolding Q' inherits the port numbering from the parent edges. 

(5) Assuming that the graph Q has a max-min LP associated with it, the max-min LP 
associated with the unfolding Q' is defined by inheritance from the parent nodes 
and edges. In particular, the type of each node (agent, constraint, objective) is the 
type of the parent node, and the coefficients associated with the edges (a^^, Ckv) are 
inherited from the parent edges. 

(6) Any two nodes of Q' with the same parent are related by an automorphism of Q' . In 
particular, any deterministic local algorithm in the port numbering model must give 
the same output on any two nodes with the same parent. Any locally computable 
feasible solution of the max-min LP associated with Q' defines a feasible solution of 
the max-min LP associated with Q, with the same utility. 

(7) Any feasible solution of the max-min LP associated with Q defines, by inheritance, 
a feasible solution of the max-min LP associated with Q' , with the same utility. 

(8) Any locally computed feasible solution of Q' with utility at least 1/a times the utility 
of any feasible solution of Q' yields an a-approximation of the optimum of Q. 

3. Local transformations 

Consider an arbitrary max-min LP associated with the graph Q. In this section we 
carry out a sequence of locally computable transformations, with the goal of arriving at a 
more structured max-min LP. The transformations are applied in the order of presentation, 
from ^3.1l to ^3.51 We describe each individual transformation in three parts: 

(1) A description of the transformation. 

(2) Mapping a solution of the transformed instance back to the original instance. 

(3) Implications to approximability. We write a;(-) for the utility of the original instance 
and for the utility of the transformed instance. 

Appendix lA.il presents the implementation of these transformations in the port numbering 
model. Figure [2] illustrates the transformations that modify the communication graph Q. 

To avoid degenerate cases, we assume that each constraint and objective is adjacent to 
at least one agent, and every agent is adjacent to at least one constraint and at least one 
objective, that is, \Vi\ > 1, \Vk\ > 1, \Ky\ > 1, and > 1 for all i G /, /c G ET, f G V. 
Indeed, isolated constraints can be deleted, isolated objectives force the optimum of ()1.2p 
to zero, non-contributing agents can be set to zero, and unconstrained agents can be set 
to +00. Furthermore, we assume that Q is connected, as we can handle each connected 
component independently. 
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Figure 2: Local transformations in ^3.11 - 

3.1. Augmenting singleton constraints — forcing \Vi\ > 2 

The transformation. For each constraint i I with \Vi\ = 1, introduce three new 
agents, s, t, and u, two new objectives, h and i, and one new constraint, j. Let v Vi 
be the original agent adjacent to i, and let k € Ky be an objective adjacent to v. Set 
Qis = ajt = Uju = 1, Chs = ces = 1, and cm = = 2 Yl^^Vk '^kw minjg/^ a,^^ 

Mapping back. Let x' be any feasible solution of the transformed instance. We map 
this back by setting = x'^ for all original agents v (zV. 

Approximation ratio. Observe that we can always set x'^ = and Xj = = 1/2 without 
decreasing the objective value uj'{x'). Thus the optima of the original and transformed 
instances coincide, and any approximation ratio is preserved. 



3.2. Reducing the degree of constraints — forcing = 2 

The transformation [3J. Replace each constraint i G I with \Vi\ > 2 by the ('^*') 
constraints 

aiuXu + ciiyXy < 1, u,veVi, u<v. (3.1) 
Mapping back. Let x' be an arbitrary feasible solution of the transformed instance. We 
map this back to a feasible solution x of the original instance by setting 

2r' 

Xy = vGV. (3.2) 

maxje/„ \ Vi\ 

To verify that x is feasible, consider an arbitrary original constraint i G /. By the 
previous transformation, we have \Vi\ > 2. Taking the sum over all the constraints (|3.1|) 
replacing i, or, if \Vi\ = 2, considering the original constraint, we have 

i:.ey,(l^^!-l)a»4<l^d(l^.|-l)/2. 

By (j3.2p we thus have 

SiiGFi '^iv^v < YlveVi '^C'ivXy/\Vi\ < 1. 

Approximation ratio. Because the objectives are unchanged in the transformation, we 
have <-i-'(x) > 2a;'(x')/A/ by linearity and ()3.2p . Furthermore, an optimal solution of the 
original instance is a feasible solution of the transformed instance. Therefore, if x' is an 
a-approximate solution, then x is a aA7/2-approximate solution. 
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3.3. Associating a unique objective with each agent — forcing \Ky \ = 1 

The transformation. For each agent v (zV with \Ky\ > 1, replace v with \Ky\ copies of v 
as fohows. Associate each copy of v with a unique objective in Ky. Replace each constraint 
adjacent to v with \Ky \ copies of the constraint, with v replaced by a unique copy in each 
constraint. The coefficients are unchanged. 

Mapping back. Let x' be a feasible solution of the transformed instance. By symmetry 
we can assume that all copies n of f have the same value without decreasing the objective 
value u;'(x'); indeed, if the values are different, just set all copies to the maximum value. 
Mapping back is done simply by identifying the copies back to the original. 

Approximation ratio. Preserved. The optima of the original and the transformed in- 
stance coincide. 

3.4. Augmenting singleton objectives — forcing \Vk\ > 2 

The transformation. For each objective k (z K with |Vfc| = 1, let v be the unique agent 
adjacent to k. Replace v with two copies, t and n, and replace each constraint adjacent to 
V with two copies of the constraint, one containing t in place of v, and the other containing 
u in place of v. Let c^t = c^u = Cfct,/2. The coefficients are otherwise unchanged. 

Mapping back. Let x' be a feasible solution of the transformed instance. By symmetry 
we can assume that the copies of v have the same value x[ = x'^ without decreasing the 
objective value a;'(x'); indeed, if the values are different, just set all copies to the maximum 
value. Mapping back is done simply by identifying the copies back to the original. 

Approximation ratio. Preserved. The optima of the original and the transformed in- 
stance coincide. 

3.5. Normalising coefficients — forcing Ckv = 1 for adjacent k & K and v € V 

For each v £ V, let k{v) be the unique objective in Ky. 

The transformation. For each v £ V, i G I, k G K, divide aiy and Ckv by Cfe(„)„. 
Mapping back. For each v £V, multiply Xy by Cfc(„)t,. 
Approximation ratio. Preserved. 

4. Local approximation algorithm 

Throughout this section we consider a max-min LP associated with a bipartite graph 
G = {V L) I U K,E) with these properties that follow from fJ5] and ^ 

(1) the graph Q is an unfolding of a finite graph, 

(2) every agent v £V is adjacent to exactly one objective and at least one constraint, 

(3) every constraint i G / is adjacent to exactly two agents, 

(4) every objective k £ K is adjacent to at least two agents, 

(5) for adjacent k £ K and v £ V, the coefficient Cky equals 1. 
It follows that Q is countably infinite. 

Recall that k{v) is the unique objective adjacent to v G V. Let N{v) = Vfc(^) \ {v} be 
the set of other agents adjacent to this objective. For a constraint i £ I and an agent v £ Vi, 
denote by n{v,i) the unique agent other than v in Vi. See Figure [U for an illustration. 

Let i? = 2, 3, ... be a fixed parameter that will determine the local horizon and the 
approximation ratio. Let r = R — 2. 
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4.1. An upper bound via alternating trees 

A walk W = (uq, ei, ui, . . . , ei,U£) in Q is alternating if (i) for all 1 < j < j' < i with 
Uj € K and Uj' G if, there exists a j < j" < j' with Uj/' G /; and (ii) for all 1 < J < j' < ^ 
with G / and Uj' G /, there exists a j < j" < j' with Uj// G -ftT. 

Let n G y be an arbitrary agent, and consider the subgraph Au of G induced by the 
nodes reachable via alternating paths starting at u that (i) traverse the constraint k{u) and 
have length at most 4r + 3; or (ii) have length at most 1. The level of a node of Au is 
its distance to k{u), with the exception of u, which we define to have level —1, and the 
constraints adjacent to u, which we define to have level —2. See Figure [1] for an illustration. 

Associate with Au a max-min LP by restriction from the max-min LP associated with Q. 

Lemma 4.1. The graph Au is a finite tree. Moreover, 

(1) every objective in Au is at level (mod 4), 

(2) every agent in Au is at level either 1 (mod 4) or 3 (mod 4), 

(3) every constraint in Au is at level 2 (mod 4), 

(4) the leaves of Au are constraints at levels —2 and 4r + 2, 

(5) for any objective k in Au and any agent v adjacent to k in Q, the agent v, occurs 
in Au and is adjacent to k in Au- 

Proof. By induction on R using the assumptions on the structure of ^. ■ 

Lemma 4.2. The optimum value of the max-min LP associated with Au is an upper bound 
on the value of any feasible solution of the max-min LP associated with Q. 

Proof. Because Au is finite, the optimum is well defined. It suffices to show that any feasible 
solution of Q is (by restriction) a feasible solution of Au- This follows from Lemma l4.ll 
the objectives in Au are identical to those in and the constraints on variables in Au are 
either identical or relaxed from Q (at leaves or non-alternating constraints). ■ 

4.2. The optimum of Au 

The tree structure of Au enables a recursive characterisation of the optimum that 
proceeds level-wise towards u. Denote by L{u,i) the set of all nodes at level i in Au. 

Let tu be the maximum value lij > such that for all d = 0, 1, . . . , r it holds that the 
values 

?;GL(n,4r + l), (4.1) 

E.gtvmO))' ^GL(n,4(r-d)-l), (4.2) 

\niv,^fu,nM4-i^))^ d>h veL{uA{r-d) + l) (4.3) 

0<d<r, V e L{u,4:{r - d) + (4.4) 

(4.5) 

Note that the maximum exists because ()4.2p and ()4.3p can be expressed using linear in- 
equalities (by introducing additional variables), and cj = is a feasible value. In a practical 



fu,v,A^) = max(0, UJ - 
satisfy the constraints 
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implementation of our algorithm, we do not need to invoke an LP solver; a simple binary 
search for an approximation of tu is sufficient. 

Lemma 4.3. Let x be any feasible solution that achieves the objective value uj in the max- 
min LP associated with Au- Then, for all d = 0,1, ... ,r it holds that 

0<x^< f:^^^^d{oj) for all v e L{u, 4(r - d) + l) ; and (4.6) 
fuvdi^) <Xv< mina"^ for all v G L{u,4{r -d)-l). (4.7) 
In particular, t„ is the optimum utility of the max-min LP associated with Aw 



Proof. By induction on d and in order of evaluation of the recursive steps ()4.2p and ()4.3p : 
see Appendix IA.2I for a full proof. ■ 

In what follows we use the shorthand notation = /^^„^rf(t«) and f~^^^ = f~^^di^u)- 

4.3. Smoothing 

For each agent v V in Q , let s^ be the minimum of the values tu over all agents u (zV 
at distance at most 4r + 2 from v in Q. For all € ^ in ^ and all d = 0, 1, . . . , r, define 

9v,o = '^^;^<\ (4-8) 
9v,d = max(0, s„ - Ew€N(v) dtj ' (4-9) 

<d = - %niv,i)9n{v,^),d-l)^ ^ ^ (^"^O) 

Lemma 4.4. For all u £ V and all d = 0,1, ... ,r it holds that 

%4 < fu,v,d for all V G L{u,A{r - d) - 1), (4.11) 

fu,v4 < 9td M all V G L{u, 4(r - d) + 1). (4.12) 

Proof. By induction on d, using the definition of ; see Appendix IA.3I for a full proof. ■ 

Lemma 4.5. For all v €z V it holds that g^^ > and < miujg/^ a~j^ . 

Proof. Let v gV he arbitrary and choose u E N{v); such a u exists because every objective 
is adjacent to at least two agents. We have g^,. > f^vr ^ by ()4.12p and ()4.4p . Similarly, 
let n = f to obtain 

Qv^r < fv,v,r < mina",^ 

by (IOT|) and (gSl). ■ 
Lemma 4.6. For all v £ V and d = 1,2, . . . ,r it holds that g^^_i < g~^ and g^^ < g^^__i. 
Proof. By induction on d; see Appendix IA.4I for a full proof. ■ 
Lemma 4.7. For all v G V and d = 0,1, . . . ,r it holds that g^^ > 0. 
Proof. By Lemma 14.51 and Lemma [ 
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K +1 up-agent 
+1 



down- agent +1 / 
Figure 3: The weights used to assign the layers. 

Finally, each agent v outputs the value 

1 



d=0 



(4.13) 



This completes the description of the algorithm. The algorithm is local, with the local 
horizon Q{R). We proceed to show the vector x is a feasible solution, and within factor 
2(1 - 1/Ak) + e' of the optimum. 



5. Analysis 

We start by partitioning the set of agents V into up-agents and down- agents such that 

(i) every constraint is adjacent to exactly one up-agent and exactly one down- agent; and 

(ii) every objective is adjacent to exactly one up-agent. 

Associate an integer layer to each node of G as follows. First, fix an arbitrary objective 
A; G -ftT to be at layer 0. Then, determine the layer of every other node u by considering 
the unique directed path connecting k to u. The layer of u is determined by taking the 
sum of the weights of the directed edges in the path, where the weights are displayed in 
Figure El Figure [T] displays an example of a partition into up- and down-agents together 
with an assignment of layers. 

Lemma 5.1. The layers of the nodes of G satisfy the following four properties: 

(1) every objective has layer (mod 4), 

(2) every down-agent has layer 1 (mod 4), 

(3) every constraint has layer 2 (mod 4), 

(4) every up-agent has layer 3 (mod 4). 

Proof. Immediate from Figure O ■ 



5.1. Shifting strategy 



Let j = 0, 1, . . . , i? — 1 be a shift parameter. For each agent v (zV, represent the layer 
of V uniquely as 4(i?c + j) + Ad + e for integers c, d, e with < d < R — 1 and e G {—1, 1}. 
Recall that r = R — 2. Associate with v the value 



Vvij) 





,9^,1 — d 



iid = R 
iid<R 
iid< R 



1; 

2 and e 
2 and e 



(5.1) 



Observe that a down- agent has e = 1 and an up-agent has e = — 1, regardless of j. 
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For an objective k and a vector z indexed by the agents, let i^fe(z) = ^^gy^ Zy. 

Lemma 5.2. The vector y{j) is a feasible solution of the max-min LP associated with Q. 
For every objective k ^ K , it holds that Wfc(y(j)) = if k is at layer Aj — 4 (mod 4i?), and 
^fc(y(j)) > min^gVfc otherwise. 

Proof From (US]), <^M, dHO]), and <^Jij; see Appendix [XS] for a full proof. ■ 
Let us now average over all values of the shift parameter j = 0,1, . . . , R — 1 to obtain 

,„ = ix;',„o) = (ip-V (5.2) 

[REd=o9ld if IS a down-agent. 

Lemma 5.3. The vector y is a feasible solution of the max-min LP associated with Q. For 
every objective k € K, it holds that uJk{y) ^ (1 — 1/-R)min^ev'fe ■5^,. 

Proof. Follows from Lemma 15.21 ■ 
5.2. Averaging 

Associate with each agent v & V a solution y^^'' defined as follows. Choose the layers 
so that V is an up-agent; this is always possible. Let y^^^ be the value of (j5.2p . 

Lemma 5.4. The vector x is a feasible solution of the max-min LP associated with Q . 

Proof. Consider an arbitrary constraint z € /; let 14 = {v,w}. Note that whenever v is an 
up-agent if is a down-agent and vice versa. Let 

z=(yT-+yT-)/2. 

By (152]) and ({^3]) we have 



Zv = ^ ^ 2 V a..,, + —y q:.,,\= Xr. 



r 1 \ 



and Zw = Xw By Lemma [5^ the solutions y^^^ and y^'" do not violate the constraint i. 
Therefore 

a^Xy + au,Xw = a^z^ + a^iZu, = {{a^yl" + a^y]^) + {a^yl^ + a^y]^))/2 < (1 + l)/2 = 1. 
We conclude that the solution x does not violate the constraint i. ■ 
Lemma 5.5. For every objective k G K , 

Proof. Consider an arbitrary objective k in Q. Note that whenever f € is an up-agent, 
each w G N{v) is a down-agent. Let 

z = — V y^^ 

By ([521) and (fi?T3]) we have 



IVfcl 

2(|y,i-i)^" 2m\-i)^^;- \v,\ 
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for all M € Vfc. By Lemma 15? 

^ ^ 2{\Vk\ - '- 2{\Vu\ - 1) RJ -^vl 



u&Vk u&Vk 



The claim follows. 



5.3. Completing the analysis 

Lemmata 14.21 and 14.31 show that for any u G y, the value t„ is an upper bound for 
the utility of any feasible solution of the max-min LP instance associated with and so 
is s^. Lemma 15.51 therefore shows that our local approximation algorithm achieves the 
approximation ratio of 2(1 — 1/A/^)(1 + 1/{R — 1)), for the special case studied in ^ 

Together with the local transformations of ^ taking into account the increase of the 
approximation in §3.21 we conclude that the max-min LP problem admits a local algorithm 
with the approximation ratio of A/(l — 1/A;^)(1 + l/(-R — 1)); the local horizon is Q{R). 
Theorem 11.21 follows by choosing a sufficiently large R. In particular, the constants A/ and 
IS.K are not required to run the algorithm. 
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Appendix A. 

A.l. Local computability of the transformations in §0] 

Even though a description of a local algorithm often involves interleaved steps of com- 
munication and computation, it should be noted that any local algorithm with local horizon 
D can always be implemented as follows: 

(1) Each node gathers full information about its radius D neighbourhood; this is the 
local view of the agent. 

(2) Each node simulates the algorithm in its local view to determine its output. 

As pointed out in fj2l we can assume that the local view is a tree. Furthermore, in our case 
it is sufficient that each agent performs these steps - constraints and objectives do not need 
to produce any output. 

Therefore we can implement each transformation presented in this section as follows 
(with a small increase of the local horizon): 

(1) Each agent gathers its local view, up to some constant distance. This is a tree. 

(2) Each agent performs the transformation in its local view. The result is a graph, 
possibly with cycles. 

(3) Each agent unfolds the graph to obtain a tree, discarding parts that are beyond its 
local horizon. 

(4) Each agent simulates the rest of the local algorithm in this tree, and applies the 
back-mapping to determine its output. 

In some of the transformations, new agent nodes are created. In ^3.H the output of the 
new agents is not needed. In §3.31 and p. 41 existing agents can simulate their copies and 
compute the back-mapping. 

The nontrivial part is to make sure that the transformations can be performed deter- 
ministically: if the local views of agents u and v partially overlap, and both agents perform 
a transformation in the common part, the common parts must be identical after the trans- 
formation. In particular, the port numbers must be identical. In the following, we show 
how to achieve this. 

§3.11 Node k is chosen to be the node in that has the smallest port number in v. 
The port number in i for {i,s} is one larger than the port number for {i,v}. Within the 
gadget defined by s,t,u, h,i,j, we can choose some fixed port numbering. 

§3.21 §3.3^ §3.41 These transformations can be summarised as follows: pick a subgraph, 
take n copies of it, and discard some edges (e.g., in §3.31 we take copies of a subgraph induced 
by an agent v €V and the adjacent constraints 1^). We can choose the port numbers of the 
transformed instance deterministically as follows. While we are taking copies of subgraphs, 
copy the port numbers as well. This may create duplicates. However, in each of these 
transformations, a set of duplicate port numbers in a node corresponds to a set of copies of 
one subgraph (e.g., in §3.31 each copy of the subgraph induced by v and 1^ creates exactly 
one duplicate port number for each agent u ^ v adjacent to any i € /„). Furthermore, we 
can impose an ordering for the copies: in §3.21 by using pairs of port number in i; in N3.3I 
by using port numbers in v; and in §3.41 bv arbitrarily ordering otherwise identical copies. 
Therefore we can re-number the ports in each node: first, order the ports by existing port 
numbers, using the ordering of the copies to break the ties; second, assign new port numbers 
1,2,... in this order. 

§3.51 The graph is not changed, and port numbering is preserved. 
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A. 2. Proof of Lemma 14.31 

To set up the base case at d = 0, observe that by the feasibiUty of x and (|4.ip we have 

< minari = /+^Q(a;) for all v G L{u,ir + 1). 

ieiv ' ' 

Next consider any < d < r and assume inductively that 

Xw < fu,w,di^) for all weL{u, 4(r - d) + l) . 

Consider an arbitrary v £ L(n, 4(r — d) — 1). Observe that for all w £ V^'(ti) it holds that 
either v = w or w L{u, 4(r — d) + 1). If we have 

w£N(v) 

then (j4.2p implies = < x^,; otherwise ()4.2p implies 

weN{v) weN{v) 

Here the first inequality follows by the induction hypothesis, and the second inequality 
follows by assumption that x achieves the objective value to; in particular, 

^ ^ Xw — -\- ^ ^ ^ UJ . 

To complete the induction, consider any 1 < d < r and assume inductively that 
fu,w,d-i(^) ^ foi^ ah w £ L(n,4(r - {d - I)) - l). 

Consider an arbitrary v € L(u,4(r — d) + 1). Observe that n{v,i) € L(ti,4(r — (d— 1)) — 1). 
Because x is feasible, we have aj^x„ + ai^n{v,i)Xn{v,i) < 1 foi^ all i S Thus, the inductive 
hypothesis and (|4.3p imply 

To conclude that is the optimum utility of the max-min LP associated with ^i^, 
consider an optimal solution x with the objective value uj* . Observe that (14. 6p and (14. 7p 
imply (14. 4p and ()4.5p . Therefore > i-i-'*. Furthermore, > ^* would contradict with the 
assumption that x is optimal. ■ 
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A. 3. Proof of Lemma 14.41 

Consider an arbitrary u (zV. For sdl w (^V at distance at most 4r + 2 from u we have, 
by the definition of Sw , 

< < tu. (A.l) 
We proceed by induction on d. To estabhsh the base case at d = 0, observe by (14. ip 
and gSD that /+^ q = 5+0 for all v £ L{u, 4r + 1). 

Next consider any < d <r and assume inductively that 

< <d for all w G L{u, 4(r - d) + 1). 

Consider an arbitrary v € L(u, 4(r — d) — 1). Observe that for all w E V^(^) either w = v 01 
w G 4(r - d) + 1). Apply (fO]) . (|A.ip . the inductive hypothesis, and to obtain 




Finally, consider any 1 < d < r and assume inductively that 

9^^4-1 ^ fu,rn,d-i for all w e L{uA{r - {d - 1)) - 1). 

Consider an arbitrary v E L(n, 4(r — d) + 1). Observe that for all i G /{, it holds that 
n{v,i) G L(ii, 4(r — (d — 1)) — 1). Apply (14.30 . the inductive hypothesis, and (I4.10p . to 
obtain 

/tt^,d = ^Y^"-iv^ ~ "'i,n{v,i)fu,n(v,i),d-l) - "^^^^^-iv^ ~ "'i,n(v,i)9n(v,i),d-l) = 9^4' 

The induction is now complete. ■ 



A. 4. Proof of Lemma 14.61 

To set up the base case at d = 1, observe first that g^Q > by (|4.9p . By (j4.10p and 
(HSI) thus 

q^^ = mina~nl— a. , ■•,q~, x „ I < mina~^ = q^n- 

^v,L IV y i,n(v,i)^n(v,i),0 J — j^^j iv au^u 

Next consider any 1 < d < r and assume inductively that g'^^ < g^^_i- Apply (|4.9p and 
the inductive hypothesis to obtain 




~ 9v,d-i- 



Finally, consider any 2 < d < r and assume inductively that g^ < g^ Apply (j4.1Up 
and the inductive hypothesis to obtain 

q~^, = min aj~,^ ( 1 — a . . ■•,q~, ■•, , ,] < minaj^nl — a. , ..q~, , „) = q~^, -, . 

^v,d IV y t,n{v,tj-^n{v,i),d—l J — tv y i,n{v,i)^n(v,i),d—2 j Jv,d—1 

This completes the induction. ■ 
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A. 5. Proof of Lemma 15.21 

Feasibility. The vector is nonnegative by (j4.9p and Lemma 14.71 Consider an arbitrary 
constraint i £ I. By Lemma fS.ll we can represent the layer of i uniquely as 4(i?c+j) + 4(i+2 
for integers c, d with < d < R — 1. Let Vi = {v, w} with v at layer 4(i?c + j) + 4(i + 1 and 
w at layer A{Rc + j) + A{d + 1) - 1. 

First consider the case d = R — 1. Note that the layer of w is actually 4(i?(c+ 1) +j) — 1, 
that is, (i = and e = — 1 for if. By (15. 1|) and Lemma 14.51 we have 

aivUvij) + ainiVwij) = aiwgw,r < 1- 
Next consider the case d = R — 2. By (j5.ip and (|4.8p . we have 

aivUvU) + aiwVwij) = aivgtfi < 1- 
Finally consider the case d < R — 2. By (j5.ip and ()4.10p . we have 

The claim follows since i was arbitrary. 

Objectives. Consider an arbitrary objective k in Q. By Lemma l5.ll we can represent 
the layer of k uniquely as 4(i?c + j) + M for integers c, d with < d < — 1. There is a 
unique up-agent in Vfc at layer 4(i?c + j) + 4(i — 1. Denote this agent by v. The other agents 
w G A^(f) are down-agents at layer 4(i?c + j) + Ad + 1. 

First consider the case d = R — 1. By (15. ip we have 

wfc(y(i)) = + J2 ^"'(•^^ = 

Then consider the case d < R — 2. By ()5.ip and (j4.9p . we have 

weN(v) weN{v) 
> - 5+ + Yl S^r-d = Sv > min s„. 

— — u^Vk 

w£N{v) w£N{v) 

The claim follows because k was arbitrary. ■ 



